Thresholds for topological codes in the presence of loss 
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Many proposals for quantum information processing are subject to detectable loss errors. In 
this paper, we show that topological error correcting codes, which protect against computational 
errors, are also extremely robust against losses. We present analytical results showing the maximum 
tolerable loss rate is 50%, which is determined by the square-lattice bond percolation threshold. 
This saturates the bound set by the no-cloning theorem. Our numerical results support this, and 
show a graceful trade-off between computational and loss errors. 



Quantum information is delicate. Qubits can be cor- 
rupted by environmental noise, dissipation and imperfect 
logic gates. Quantum error correction [Tl[2] and fault tol- 
erant quantum computation (FTQC) [3| address these 
problems, enabling quantum information to be stored 
and manipulated in spite of physical errors. 

Errors may be classified as computational errors, in 
which the state of the qubits remains within the compu- 
tational basis, or as losses, in which physical qubits (e.g. 
a photon) are lost from the computer. More generally, 
any detectable leakage process taking a qubit out of the 
computational basis can be treated as a loss error. 

Losses are both both detectable and locatable, suggest- 
ing they should be easier to rectify than computational 
errors. Indeed, quantum communication channels can tol- 
erate a higher rate of loss {pioss < 0.5) than depolarisa- 
tion {pdepoi < 1/3) [4]. Furthermore, a quantum compu- 
tation scheme has been developed which works as long 
as pioss < 0-5, saturating this bound [5 . This thresh- 
old is much less restrictive than the largest estimates for 
the computational error threshold, Pcom ^ 10~^ [3 
however it is not clear how the scheme performs in the 
presence of both loss and computational errors. Dawson 
et al. ^ have considered an error model which contains 
both loss and computational errors, finding that FTQC 
is possible provided pioss ^ 3 x 10"^ and Pcom ^ 10""^. 
It is natural to ask whether alternative FTQC schemes 
can be found, which are tolerant to both types of error, 
and have less restrictive thresholds. 

In this paper we partially address this question by con- 
sidering the effect of both qubit losses and computational 
errors on Kitaev's surface codes [9 . Surface code quan- 
tum memories are robust against computational errors, 
with a threshold of Pcom < 0.104, and are an ingredi- 
ent in an FTQC scheme with a high computational error 
threshold, Pcom ^ 10~^ [7]. We show that surface codes 
can be made robust against both computational and loss 
errors. Furthermore, we show that the threshold for the 
loss rate on the surface code is 0.5 (when Pcom = 0) satu- 



^TMS and SDB contributed equally to this work. 
* stace@physics.uq.edu.au 
seandbarrett@gmail.com| 




FIG. 1: (a) Physical qubits (arrows) reside on the edges of a 
square lattice (dashed). Also depicted are a plaquette opera- 
tor, a star operator, the logical Z operator and the logical X 
operator, (b) In the event of a qubit loss, an equivalent logical 
operator Z can be routed around the loss. 



rating known bounds. Here, we focus on quantum memo- 
ries with idealised, perfect measurements. We also discuss 
how the insights gained in this work generalise to FTQC 
[7 which should therefore be robust to loss, computa- 
tional errors and imperfect measurement. This is relevant 
to various implementations of cluster-state computation, 
including optical and atomic ensembles |lTl . 

For the purposes of analysis, the error model we con- 
sider is local and uncorrelated. Each physical qubit is 
lost with probability pioss • Losses are presumed to be de- 
tectable: a projector onto the computational basis of a 
given qubit, 11^ = |0)i(0| + |l)i(l|, is an observable in- 
dicating only whether the state of the qubit has leaked 
out of the computational basis. The remaining qubits are 
subject to independent bit-flip {X) and phase (Z) errors, 
each with probability Pcom- Both errors are handled in an 
analogous way in the surface code, so here we confine our 
attention to X errors, noting that the thresholds for Z er- 
rors will be identical. Aside from these errors, we assume 
other quantum operations (e.g. syndrome measurements) 
can be implemented perfectly. 

Kitaev's surface codes are defined by a set of stabilisers 
acting on a set of physical qubits that reside on the edges 
of a square lattice [13] . The stabiliser group is generated 
by plaquette operators, which are products of Z opera- 
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FIG. 2: (a) Lattice with two lost qubits (crosses). Represen- 
tative qubits and plaquettes are labeled 1 to 8, and A to E 
respectively, (b) Plaquettes sharing a lost qubit ({A, B} and 
{C, D}) become superplaquettes AB and CD, and may be mul- 
tiply connected (i.e. share more than one qubit). (c) Degraded 
lattice showing superedges (thick lines), (d) Restored lattice 
with zero weight edges (dotted) and irregular weights (thick). 



tors acting on qubits within a plaquette, Pp = 
and star operators, which are products of X operators 
acting on qubits within a star, Ss = (8)jGs depicted 
in Fig.[T]^. The stabilisers commute, and the code space 
{|C)} is a simultaneous +1 eigenstate of all stabilisers. 
If the L X L lattice has periodic boundary conditions, 
there are 21/^ physical qubits and 2(L^ — 1) independent 
stabilisers. The two remaining degrees of freedom are ca- 
pable of encoding two logical qubits, Qi (i G {1,2}). A 
logical Zi {Xi) operator corresponds to a product of Z 
{X) operators along a homologically nontrivial cycle (i.e. 
spanning the lattice), shown in Fig. [l] Zi and Xi com- 
mute with the stabilisers, but are not contained within 
the stabiliser group. 

The logical operators are not unique; e.g. multiply- 
ing Zi by a plaquette stabiliser yields a new oper- 
ator Zi = PpZi that acts identically on the code 
space, Zi\C) = PpZi\C) = ZiPp\C) =Zi\C). Thus there 
are many homologically equivalent cycles with which 
to measure each logical qubit operator, as shown in 
Fig. [ijb). This redundancy allows us to obtain the loss 
threshold for the case Pcom = 0: if only a few physical 
qubits are lost, it is likely that each logical operator can 
be reliably measured by finding a homologically nontrival 
cycle that avoids all lost qubits. 

If pioss is too high, there is likely to be a percolated 
region of losses spanning the entire lattice, in which case 
there are no homologically nontrivial cycles with which 
to measure the logical operators. As L ^ oo, there is a 
sharp boundary between recoverable and non-recoverable 
errors corresponding to the bond percolation threshold on 
the square lattice [14 : ioi pioss < 0.5 loss recovery almost 



surely succeeds, whereas for pioss > 0.5 loss recovery al- 
most surely fails. Notably, this threshold saturates the 
bound on pioss imposed by the no-cloning theorem [4]. 

The case pioss = and Pcom > has been well studied 
|9l [m [16] . Briefly, physical bit-flip errors lead to logical 
bit-flip (Xi) errors but not logical phase errors, and vice- 
versa. An error chain^ is a set of lattice edges (i.e. 
physical qubits) where a bit-flip error has occurred. The 
plaquette operator eigenvalues change to —1 only at the 
boundary^ dE^ of the chain. Measuring the plaquette op- 
erators therefore yields information about the endpoints 
of connected sub-chains of ^. If £^ crosses Zi an odd num- 
ber of times, then the logical qubit suffers a Xi error. 
These errors may be corrected if the endpoints, dE^ can 
be matched by a correction chain, E' ^ such that the closed 
chain C = E ^ E' crosses Zi an even number of times, 
i.e. C is homologially trivial. The error rate below which 
the correction chain E' may be successfully constructed 
is closely related to the phase boundary of the random- 
bond Ising model (RBIM) [151 [IT]. If Pcom < PcQ = 0.104 
[18^, then in the limit L ^ oo, the most probable chain, 
Cmax = E -\- E'^^^ , is almost surely homologically trivial 
and recovery succeeds. If Pcom > PcO^ then in the limit 
L ^ oo, the chain is homologically trival only 25% of the 
time, and recovery fails. 

We can think of the above results as endpoints of a 
'boundary of correctability': (piossiPcom) = (0.5,0) and 
(0,0.104), respectively. In what follows, we demonstrate 
that toric codes (and planar codes, by extension) are ro- 
bust against both loss and computational errors, with a 
graceful tradeoff between the two. We first describe how 
losses can be corrected by forming new stabiliser genera- 
tors, which are aggregations of plaquettes or stars, called 
superplaquettes and superstars, respectively. The super- 
star and superplaquette eigenvalues then reveal the error 
syndromes, and a perfect matching algorithm is used to 
find an error correction chain E' . We illustrate the ef- 
ficacy of the single round error correction protocol by 
calculating numerically the boundary of correctability in 
the {pioss.Pcom) parameter space. 

Consider the lattice shown in Fig. |2|a) which is dam- 
aged by the loss of two physical qubits, marked by the 
crosses. The loss of qubit 3 affects two plaquette sta- 
bilisers: Pa = ZiZ2ZsZ^ and Pg = ZsZ^ZqZ^, ren- 
dering them unmeasurable. However, the superplaquette 
Pab = ^A^B = ZxZ^Z^Zr^Z^Z^ is independent of the 
qubit at site 3, so stablises the remaining qubits. Without 
errors, Pab has an eigenvlaue of +1. An error chain end- 
ing within the superplaquette AB changes the eigenvalue 
of Pab to —1. It follows that the syndrome associated 
with a superplaquette is determined by the parity of the 
number of error chains that cross its boundary. The fact 
that superplaquette operators yield syndrome informa- 
tion with which to construct an error correction chain, E' ^ 
is the basis for our loss-tolerant error-correction scheme. 

In general, given any set of lost qubits, we can form a 
complete set of stabilisers on the damaged lattice in the 
following way: for each lost qubit which participates in 
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neighbouring (super)plaquettes Pq and P^, we form the 
superplaquette operator PqPq, which is independent of 
Zq. In the same way, we form superstar operators from 
products of star operators. As discussed earlier, we can 
also form new logical Xi and Zi operators by deforming 
the original logical operators to conform to the bound- 
aries of newly formed superplaquettes. 

We note that in Fig. [2|a), there is a damaged pla- 
quette operator Zj = Z^Pa = Z1Z2Z4 (or, equiva- 
lent ly ZsPs = Z^ZqZj) associated with the lost qubit 
3, which commutes with all the newly formed stabiliser 
generators on the damaged lattice, but whose eigenvalue, 
±1, is indeterminate. Likewise, the damaged star oper- 
ator Xj = X4X7X8 also has indeterminate eigenvalue 
and commutes with the new stabilisers on the damaged 
lattice. Having indeterminate eigenvalues, Zj and Xj^ 
which mutually anticommute, define a two-dimensional 
degree of freedom in an uncertain state. They therefore 
describe a completely mixed junk qubit, J'^ which is a 
consequence of the entanglement between the lost qubit 
and the remaining qubits [19^. Since Zj and Xj each 
commute with the new stabilisers, and with the deformed 
logical operators, the junk qubit is in a product state with 
the logical qubits: \i^){ip\qi \^){4'\q2 ^^^J^/^, and so the 
loss does not affect the logical qubits. 

When analysing the pattern of syndromes on the pla- 
quettes and superplaquettes, we construct a new graph, 
depicted in Fig. |2|b), in which a (super)plaquette is rep- 
resented by a node, and (super)plaquettes share a bond 
on the new graph wherever they share a physical qubit 
in common. Thus Pab and Pqd share the qubits 2 and 5, 
and this is represented as two edges between the super- 
plaquette nodes labelled AB and CD. 

The error correction syndrome, dE^ arising from an 
error chain on the graph in Fig. [2|b) is determined by 
the (super)plaquettes that have an eigenvalue of —1. To 
correct the errors, we follow the procedure described by 
Harrington et al [151 ES] to find the most likely error 
chain giving rise to dE. The probability of a given error 
chain is modified by the presence of losses. With no loss, 
the probability of an error on a qubit, i = {P, P'}, be- 
tween two neighbouring plaquettes P and P', is uniform, 
P£ = Pcom- With loss, superplaquettes may share multiple 
physical qubits in common, as shown in Fig. |2|b), where 
superplaquettes AB and CD have qubits 2 and 5 in com- 
mon. A non-trivial syndrome arises only if either qubit 2 
or qubit 5 suffers an error, but not both. By extension, 
for a pair of neighbouring superplaquettes, £ = {P, P'}, 
sharing physical qubits, a non-trivial syndrome arises 
only if there are an odd number of errors on the ni qubits, 
which happens with probability 

m odd 

We therefore degrade the graph shown in Fig. |2|b) , re- 
placing multi-edges (i.e. several shared physical qubits) 
whose error probabilities are uniform, with single su- 
peredges whose error rates depend on the number of 



physical qubits shared between neighbouring superpla- 
quettes. This degraded lattice is shown in Fig. |2|c), in 
which there are no multi-edges, but the error probabili- 
ties are no longer constant. 

On this degraded lattice, we may now assign a proba- 
bility for any hypothetical chain E' = E -\- where C is 
a closed chain. This probability, which is conditioned on 
the measured syndrome, dE is |16] 

P{E'\dE) = AfHe^'-^' , where = | ;1 ^ ^ ^ 

for a chain N = Y\\fi \/p^(l — Pi) is a normalisation 
constant and e^^^ = — 1. 

The chain E' that maximises P{E'\dE) also minimises 
^i^E' "^^^^ minimisation may be accomplished us- 
ing Edmonds' minimum- weight, perfect-matching algo- 
rithm [20]. For pioss = this simply minimises the total 
metropolis length of the matching path, and is the same 
procedure implemented in previous studies [151 US] • For 
Pioss > the edge weights are not uniform, since p£ 
depends on the number of physical qubits, n^, shared be- 
tween adjacent superplaquettes. 

For the purposes of simulation, it is easier to deter- 
mine homology classes on a square lattice, rather than 
the degraded lattice, exemplified in Fig. [2|c). We there- 
fore restore the square lattice by dividing superplaquettes 
into their constituent plaquettes in the following way: 
(1) an edge between two plaquettes within a single su- 
perplaquette is assigned a weight of zero, (2) an edge 
between plaquettes in two neighbouring superplaquettes 
is given the weight of the superedge in the degraded lat- 
tice, as illustrated in Fig.|2|d). These transformations do 
not change the weighted distance between any pair of 
syndromes, and so a minimum- weight perfect matching 
on the restored lattice is also a minimum-weight perfect 
matching on the degraded lattice. Determining the ho- 
mology class is then accomplished by counting crossings 
of vertical and horizontal test lines in the dual lattice. 

In order to test the efficacy of our loss-tolerant error 
correction scheme, we generate random losses on a peri- 
odic lattice with rate pioss- On the remaining qubits we 
generate a set of errors, E^ with rate Pcom- Applying Ed- 
monds' algorithm to dE on the weighted lattice yields the 
maximum-likelihood error correction chain, E' . The ho- 
mology class of the chain E^E' then determines whether 
error correction was successful. 

For each value of pioss we simulate the protocol for 
different values of Pcom on lattice sizes, L = 16, 24 and 
32. For given values of Pcom and L, the failure rate, Pfaih 
is calculated by averaging over 10^ trials. Following [15], 
we seek a threshold, pt (depending on pioss)^ such that 
dpfaii/dL < when Pcom < Pt, and dpfaii/dL > when 
Pcom > Pt- That is, for each value of pioss, we fit the 
simulated failure rate to a universal scaling law Pfaii = 
f[x] a-\- 6x, for small x), where x = {pcom — Pt)L^^^^ -> 
with fitting parameters pt^i'o^CL and b. 

Fig. [3] is the central result in this paper, and shows 
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FIG. 3: Correctability phase diagram. The shaded region is 
correctable in the Umit L ^ oo. The threshold, pt, is calcu- 
lated by fitting the universal scaling law Pfaii = f[{pcom — 
pt)L^/^°]. The curve is a quadratic fit to the points for which 
Ploss < 0.4 (where universal scaling is unaffected by the finite 
lattice size). It extrapolates through (pioss,Pt) — (0.5,0). 



Pt as a function of pioss- For the points pioss ^ 0-4, the 
universal scaling law provides a good fit to the simulated 
results, so pt is well defined, and ~ 1.5 consistent 
with the scaling exponent of the RBIM universality class 
[15]. A quadratic fit through these points yields a curve 
that extrapolates through (0.5,0). This curve represents 
the boundary of correctability: if (pioss^Pcom) is in the 
shaded region then the failure rate decreases to zero as 
L increases. Importantly, this boundary passes through 
the known bounds at pioss = and 0.5, demonstrating 
that the protocol is very robust against loss. 

For Ploss ^ 0.425, the universal scaling assumption 
breaks down (best-fits for pt are still shown), and the 
points in Fig.jsjlie below the quadratic extrapolation (but 
still attaining the point (0.5, 0)). This is attributed to the 



fact that for pioss ^ 0.425, the largest superplaquette on 
an L < 32 X 32 lattice occupies approximately half of the 
lattice sites [21 , so finite size effects dominate. 

The protocol described in this paper for dealing with 
losses in a surface code relies on several important prop- 
erties of the stabilisers. Firstly, if a physical qubit q in 
the logical qubit operator chain Zi is lost, then there is 
a plaquette P^, such that Zi = ZiPq is independent of 
q (likewise Xi Xi). Thus, logical operator chains can 
be rerouted around the lost site. Secondly, there is an- 
other plaquette such that the superplaquette PqPq 
is independent of q. Thus, superplaquettes may be con- 
structed to locate the endpoints of error chains (likewise 
for superstars). Thirdly, newly formed junk qubits are 
uncorrelated with the logical qubits. 

These properties are satisfied by a number of re- 
lated protocols, including the fault-tolerant planar code 
[I5t [16] , long-range entanglement creation in a noisy sys- 
tem p2], and in Raussendorf's topological, fault-tolerant, 
cluster-state quantum computation scheme [7 , which 
borrows a number of ideas from anyonic quantum compu- 
tation. In these protocols, the syndrome operators exhibit 
simple cubic symmetry, whose bond percolation thresh- 
old is psc ^ 0.249 [14 . We therefore expect that the 
region of correctability for these schemes includes the 
points {pioss.Pcom) = (0,0.029) and (0.249,0). 

We have demonstrated that surface codes are robust 
to noise arising from both errors and losses. The cor- 
rectable phase in (pioss^Pcom) space includes the known 
results, and we have shown that for a model of uncor- 
related noise, there is a significant fraction of parameter 
space in which the surface code is robust to both loss 
and computational errors. The approach described here is 
applicable to other systems, including full, fault-tolerant 
quantum computation, as well as correlated noise models. 
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